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Abstract 

We present a new five-parameter class of Ricci-flat solutions in four dimensions with Eu¬ 
clidean signature. The solution is asymptotically locally fiat (ALE), and contains a finite 
asymptotic NUT charge. When this charge is sent to infinity, the solution becomes asymptot¬ 
ically locally Euclidean (ALE), and one in fact obtains the Ricci-flat Plebahski-Demiahski 
solution. The solution we have found can thus be regarded as an ALE generalisation of 
the latter solution. We also show that it can be interpreted as a system consisting of two 
touching Kerr-NUTs: the south pole of one Kerr-NUT touches the north pole of the other. 
The total NUT charge of such a system is then identified with the asymptotic NUT charge. 
Setting the asymptotic NUT charge to zero gives a four-parameter asymptotically fiat (AF) 
solution, and contained within this subclass is the completely regular two-parameter AF 
instanton previously discovered by the present authors. Various other limits are also dis¬ 
cussed, including that of the triple-collinearly-centered Gibbons-Hawking solution, and an 
ALE generalisation of the C-metric. 


1. Introduction 


Exact solutions have played an important role in the development of Einstein’s general 
theory of relativity, describing important predictions such as black holes. The hrst known 
exact solution was discovered by Schwarzschild in 1918, although it was only realised much 
later that it describes a static black hole. The rotating generalisation of the Schwarzschild 
solution was discovered by Kerr in 1963. It is a two-parameter class of solutions, describing 
a rotating mass in general relativity. 

The Schwarzschild solution has another well-known generalisation, the so-called Taub- 
NUT solution [T]. It possesses some peculiar properties, such as the presence of a so-called 
Misner string. Indeed, the Taub-NUT solution was used by Misner as “a counterexample to 
almost anything” [2]. When a rotational parameter is added, it generalises to the Kerr-NUT 
solution, which is contained within the so-called Carter-Plebahski solution ill]. The latter 
solution is more general, containing electric and magnetic charges, as well as a cosmological 
constant. But in the present paper, we will only be interested in its Ricci-flat limit, i.e., the 
Kerr-NUT solution. 

A black hole can also accelerate in general relativity. The solution that describes such 
a black hole is known as the C-metric. The C-metric was found quite early on, but its inter¬ 
pretation as an accelerating black hole was only made clear after the work of Kinnersley and 
Walker [5] in 1970. Kinnersley also found its generalisation with rotation and NUT charge 
|6]. In 1976, Plebahski and Demiahski obtained the most general black hole solution, and 
presented it in a remarkably compact form |7| . It contains seven parameters: mass, rotation, 
NUT charge, acceleration, cosmological constant, electric and magnetic charges. Restrict¬ 
ing to the Ricci-flat class, the cosmological constant and electric and magnetic charges are 
removed. We shall call such a class of solutions the Ricci-flat Plebahski-Demiahski solu¬ 
tion, although it is perhaps more appropriate to call it the Kinnersley solution. It is a 
four-parameter class of solutions. 

In this paper, we focus on solutions with Euclidean signature, i.e., those with all¬ 
plus signature. Interestingly, all the above-mentioned black-hole solutions have Euclidean 
sections. Many calculations in the black-hole space-times, in particular those involving 
quantum helds, are hrst done in their Euclidean sections and then analytically continued 
back to Lorentzian signature. Euclidean solutions are of interest in their own right as well. 
In the Euclidean path integral approach to quantum gravity [8], all the Euclidean metrics 
on a given manifold with hxed boundary conditions are integrated over. Those metrics 
satisfying the Einstein equations are thus the stationary phase points of the path integral. 
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We should, however, point out that it is not guaranteed that a Euclidean solution can be 
obtained from a Lorentzian one by Wick rotation, or vice versa. For example, the self-dual 
class of Euclidean solutions known as the multi-Taub-NUT solution in the Gibbons-Hawking 
ansatz [9i[in] has no Lorentzian section. Nevertheless, as mentioned, large classes of solutions 
do have Euclidean sections, and these will be our main focus from now on. From this point, 
we will implicitly refer to its Euclidean section when we talk about a solution. For example, 
when we refer to the Kerr-NUT solution, we really mean the Euclidean Kerr-NUT solution. 

In the Euclidean section, the time-translational symmetry becomes either a translational 
or a rotational symmetry. The Euclidean spaces now possess two commuting Killing vectors, 
one corresponding to the Euclidean time flow, and the other to the usual rotation. The black- 
hole horizons, as well as the so-called acceleration horizons when the acceleration parameter 
is present, become axes of the space. This means that some linear combination of the two 
Killing vector helds vanishes at each of these horizons or axes. The classihcation of Euclidean 
solutions based on the hxed-point sets of the Killing vector helds was carried out by Gibbons 
and Hawking and more recently in terms of the so-called rod structure by the present 
authors [12]. 

In the rod-structure formalism, the various axes of the space are known as rods. Each 
rod has a direction, which is dehned to be the normalised Killing vector held which vanishes 
along that rod. The points at which adjacent rods meet are called turning points. The reader 
is referred to [12], and references therein, for more details of the rod-structure formalism, 
and for explicit examples of rod structures of a number of Euclidean solutions. 

Euclidean solutions can also be classihed by their behaviour at asymptotic inhnity as 
being asymptotically hat (AF), asymptotically locally hat (ALE), asymptotically Euclidean 
(AE) or asymptotically locally Euclidean (ALE) [13]. A solution is said to be AF, if at 
inhnity the metric approaches the form: 

ds^p = dr^ -I- dr^ -|- r‘^{d6‘^ + sin^ 6 d0^), (1.1) 

sufficiently fast. Solutions in this class include Schwarzschild and Kerr. On the other hand, 
if the metric approaches 

ds^LF = + 2n cos 6 d0)^ -|- dr^ -|- r^(d6*^ -|- sin^ 6 d0^), (1.2) 

sufficiently fast at inhnity, it is said to be ALE. Here the so-called NUT charge n manifests 
itself in the asymptotic behaviour of this class of metrics. Solutions of this type include 


3 


ALF metric 

ALE limit 

Number of 
turning points 

Taub-NUT 

hat space 

1 

Kerr-NUT 

double-centered 

Gibbons-Hawking 

2 

n-centered 

Taub-NUT 

n-centered 

Gibbons-Hawking 

n 


Table 1: Known examples of ALE metrics as special limits of ALF metrics with infinite NUT 
charge. The last column refers to the number of turning points in the rod structures of the 
corresponding ALF and ALE metrics. 


Taub-NUT and Kerr-NUT. If the metric at inhnity instead approaches 

ds^LE = cos^ 9 dr^ + dr^ + r^{d6‘^ + sin^ 6 d0^), (1.3) 

sufficiently fast, it is called AE or ALE, depending on whether the asymptotic constant r- 
surface is identihed. We refer to these two classes as ALE collectively when global geometry 
is of no concern; in fact, AE is just a special case of ALE with a trivial identihcation group. 
Flat space is of course AE. A well-known example of an ALE space is the double-centered 
Gibbons-Hawking space, with the Eguchi-Hanson space as a special case. 

If one is only interested in local geometry, AE metrics are a special limit of ALF metrics 
with vanishing NUT charge. This is clear from the asymptotics fll.ip and fll.2p . On the other 
hand, AE and ALE metrics are a special limit of ALF metrics with infinite NUT charge. 
The latter correspondence is not manifest in the asymptotics fll.2p and fll.3p alone, but it is 
true for all explicitly known solutions (c.f. Table 1). For example, the inhnite NUT charge 
limit of Taub-NUT gives flat space. In this sense, ALF metrics are the most general class of 
solutions. 

It turns out that the C-metric is ALE: the acceleration horizon now becomes one of the 
asymptotic axes in the Euclidean section. The more general Ricci-flat Plebahski-Demiahski 
solution is also ALE, as can be checked explicitly. One is then naturally led to the question 
of whether the latter metric admits an ALF generalisation. Such a generalisation, if it 
exists, would contain an asymptotic NUT charge n, such that when n is taken to inhnity, 
the Ricci-hat Plebahski-Demiahski solution is recovered. 

The aim of this paper is to explicitly present this ALF generalisation of the Ricci-hat 
Plebahski-Demiahski solution, and to analyse its properties and various limits. Like the 
Ricci-hat Plebahski-Demiahski solution, this new solution has three turning points in its 
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rod structure. It can be regarded as the non-self-dual counterpart of the triple-collinearly- 
centered Taub-NUT solution, just as the Kerr-NUT solution is the non-self-dual counterpart 
of the double-centered Taub-NUT solution in the two-turning-point case. In analogy with 
the Kerr-NUT solution, we will see that this new solution in fact contains the self-dual 
triple-collinearly-centered Gibbons-Hawking solution as another ALE limit. 

It should be emphasised that the asymptotic NUT charge n present in the new ALF 
solution that we have found, is different from the NUT-charge parameter that the Plebahski- 
Demiahski solution is known to possess. In particular, in the limit when n is taken to 
inhnity, the latter parameter is still present in the solution and remains hnite. The difference 
between these two types of NUT charges can be understood in terms of the rod structure 
of the solution. It turns out that the asymptotic NUT charge is a property relating the two 
asymptotic rods (the hrst and fourth rods of the rod structure), while the traditional NUT 
charge is an analogous property relating the hrst and third rodsj^ 

The paper is organised as follows. We begin in Sec. 2 by presenting the general metric 
together with an analysis of its symmetries. We also briehy discuss the construction of 
the solution. The main geometrical properties of the solution are analysed in Sec. 3. This 
includes the conditions necessary to ensure the correct metric signature and absence of 
curvature singularities, as well as its asymptotic structure and rod structure. In Sec. 4, we 
show how four important limits of the solution can be obtained. They all share the property 
that the number of turning points in the rod structure remains hxed at three. In Sec. 5, we 
provide an alternative interpretation of the solution as a system consisting of two touching 
Kerr-NUTs, while in Sec. 6, we discuss the interpretation of this solution in Kaluza-Klein 
theory as a system of three collinearly centered monopoles. The paper ends with a discussion 
of possible generalisations and future problems. There is also an appendix, which contains 
details of how the two- and one-turning-point limits of the solution are taken. 


2. The metric 

The solution we have found can be compactly written in the following form: 

2 {FdT + Gdct>f kH fdx^ dy^ d A 

^ {x - y)HF {x-yY\X Y kF ^ 

^The NUT-charge parameter in the Plebahski-Demiahski solution is so called because it reduces to the 
NUT-charge parameter in the Kerr-NUT solution when the zero-acceleration limit is taken. In this limit, 
the fourth rod disappears from the rod structure, and the first and third rods become the asymptotic rods. 
Thus, it is only in this limit that the NUT-charge parameter in the Plebahski-Demiahski solution has its 
usual interpretation as in 
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H = {ux + y) [{vx - y) (ai - a 2 ,xy) - 2(1 - z/) (oq - a^x^y "^)], F = y^X - x^Y , 

G = {u'^ao + 2ua3y^ + 2ua4y^ — a4y^)X + (oq — 2uao — 2uaix — i'‘^aix'^)Y , 

X = uq + aix + a2x‘^ + asx^ + 040 ;^, Y = uq + aiy + a2y‘^ + a^y^ + 04 ?/^, (2.1) 

in the coordinate system {T,(j),x,y). There are seven apparent parameters in this solution: 
an overall scale factor fc, a parameter u, and hve arbitrary coefficients ao ,,,,,4 for the quartic 
polynomial X (or Y). 

2.1. Symmetries 

It can be checked that the solution fl2.1l) has the following two continuous symmetries: 
• Scaling symmetry: 


Oj —)■ Oj/c*, X —)■ cx , y ^ cy , 


( 2 . 2 ) 


• Parameter symmetry: 


flj —^ CQii , <p —^ <p! C 


for an arbitrary non-zero constant c, and i = 0,..., 4. 


(2.3) 


This implies that two out of the five coefficients ao ,...,4 are actually redundant. So fl2.ip is in 
fact a five-parameter solution. There are also two discrete symmetries present: 


• Swapping symmetry: 


z/—>-l/z/, X-H-y , (r, 0) —)■ i(r/z/— 2020, —z/(/)), fc —)■ — fcz/^, (2.4) 

where i is the imaginary unit; 

• Inversion-swapping symmetry: 

11 , , 
X —y — , y —y — , og t-y 04 , oi i-y 03 . (2-5) 

y X 

One can use these symmetries to narrow the ranges of the coordinates and parameters 
of the solution. For example, the swapping symmetry can be used to restrict the range of z/ 
to —1 < z/ < 1 without any loss of generality. This will be used in our discussion below. 
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2.2. Alternative parameterisation 

It is often convenient to reparameterise the solution in terms of the roots of the poly¬ 
nomial X as 

X = — xi){x — X 2 ){x — X 3 ){x — Xi) , (2.6) 

by dehning 

Oo = aiXiX 2 X^Xi , 02 = ai{XiX 2 + XiX^, + XiX^ + X 2 X 3 + X 2 X 4 ^ + X 3 X 4 ) , 

ai =-ai{xiX 2 X 3 +xiX 2 Xi +X 1 X 5 X 4 +X 2 X 3 X 4 ,), =-a 4 {xi + X 2 + xs + X 4 ). (2.7) 

Then the scaling symmetry implies that the following substitutions: 

04 —)■ 04 / 0 ^ , Xi —)■ cXi , X ^ cx , y ^ cy , 0 —)■ ( 2 . 8 ) 

is a symmetry, so the relevant parameters are the relative ratios, rather than the dehnite 
values, of the roots of X. The parameter symmetry implies that the quartic coefficient 04 of 
X can be set to an arbitrary non-zero value. 

2.3. ISM construction 

The solution was constructed using the inverse-scattering method (ISM) [HI |T5]. It 
was originally obtained by applying a three-soliton transformation on the Euclidean triple- 
collinearly-centered Schwarzschild solution. The three so-called Belinski-Zakharov (BZ) pa¬ 
rameters were then hxed to eliminate the three corresponding turning points where the 
soliton-transformation was performed. One possible choice of these BZ parameters then led 
to the above solution (the other choice leads to the triple-collinearly-centered Taub-NUT so¬ 
lution). It was subsequently realised that this solution can also be generated by applying a 
three-soliton transformation on the Euclidean double-Schwarzschild solution, and then elim¬ 
inating one turning point by hxing the corresponding BZ parameter. The latter construction 
was in fact previously carried out in [16], from which the new AF gravitational instanton 
was extracted as a special case. 

The remaining major challenge was to then cast the solution in a compact form. We 
eventually succeeded by using the C-metric-like coordinates {x, y). The basic idea was to ab¬ 
sorb as many parameters as possible into the structure function X. A Mobius transformation 
can be performed on (a:, ?/), and using this freedom we hxed one factor of H as vx + y. 

The details of the above ISM constructions and simplihcations will be presented else¬ 
where. 
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3. Analysis of the geometry 


3.1. Ranges of coordinates and parameters 

To simplify our analysis, we need to narrow the ranges and remove certain redundancies 
of the parameters. We use the swapping symmetry to restrict the range of u, and choose a 
negative X 2 using the scaling symmetry: 

— I < n < 1 , X2 < 0 . (3.1) 

These ranges shall be assumed throughout the paper. X 2 can be further gauge-hxed to be — 1 
due to the scaling symmetry, and can be hxed to be 1 due to the parameter symmetry: 

X2 = -I, 04 = 1 . (3.2) 

We shall use these gauge-hxed values whenever it is convenient. 

We now deduce the appropriate ranges of the coordinates {x,y). To ensure that the 
signature of the space does not change in the region of interest, we note that x must lie 
between a pair of adjacent roots of X, and similarly for y. Furthermore, these two pairs 
of adjacent roots must be such that X has opposite sign to F, to ensure a positive- or 
negative-dehnite signature. In particular, this implies that x and y lie in different ranges. 

We also require that the region of interest is non-compact with an asymptotic infinity. 
As can be seen from the metric fl 2 .ip . asymptotic infinity is reached when x = y. Given 
that X and y necessarily lie in different ranges, there must exist a common boundary at 
which they coincide. Without loss of generality, we take this common boundary to be at 
X = y = X 2 . We further assume that x ranges between the two roots X 2 and X 3 , and that y 
ranges between the two roots xi and X 2 . 

At this stage, we note that the region of interest can be visualised as a rectangle or 
“box” in a two-dimensional plot with x and y as the axes. This box will have four sides, 
corresponding to the boundaries of the ranges of x and y. Asymptotic infinity itself is 
represented by the diagonal line x = y in this plot. The box will then have to touch the 
point X = y = X 2 on this line at either its lower-right or upper-left corner. 

Depending on the values of Xi and X 3 relative to X 2 , the box may either be finite in 
extent, or up to two of its sides may extend to infinity. For example, consider a box touching 
the point x = y = X2 nt its lower-right corner. If we have the ordering x^ < X2 < xi, the box 



will be finite in extent. On the other hand, if X 2 < the left side of the box will extend 
to infinity and “wrap around” to the right side of the plot, ending aX x = x^. Similarly, if 
xi < X 2 , the top side of the box will extend to infinity and wrap around to the bottom side of 
the plot, ending at y = Xi. We thus have four qualitatively different cases: (i) x^ < X 2 < Xi; 
(ii) X 2 < X 3 and X 2 < Xi; (iii) X 3 < X 2 and Xi < x 2 ] and (iv) xi < X 2 < X 3 . 

We next turn to the constraints set by requiring that the region of interest is free of 
curvature singularities. We note that there are possible curvature singularities whenever H 
vanishes. As can be seen from the expression of H in fl 2 .ip . this occurs along the straight 
line y = —ux, as well as on non-trivial curves where the second factor of H vanishes. We 
also note that there are curvature singularities at the four points {x = 0 ,y = ±cxd) and 
(x = ±cxo, j/ = 0). 

We now show that a box touching the point x = y = X 2 at its lower-right corner will 
always contain a curvature singularity. As explained above, there are four different cases 
to consider, labelled by (i)-(iv). It turns out that cases (iii) and (iv) can immediately be 
eliminated, since we must have X 2 < xi < 1 in order for the box not to touch the line 
singularity y = —vx for —1 < z/ < 1. Turning to case (ii), note that we require xi < 0 in 
order for the box not to contain the singularity at (x = — 00 , y = tS). If this holds, then case 
(ii) is in fact equivalent to case (i) under the inversion-swapping symmetry. Thus we are 
left with case (i) to consider. It can be checked that in this case, H will necessarily change 
sign somewhere in the box. Indeed, if one demands that the box avoids the line singularity 
y = —z/x, it will inevitably contain part of the non-trivial singularity curves. 

Hence we only need to consider a box touching the point x = y = X 2 at its upper-left 
corner. The analysis of this case can be carried out in a similar fashion. Since {x = 0,y = 
— 00 ) is a curvature singularity, any box which does not contain this point has xi < X 2 , or can 
be transformed to such a form by applying the inversion-swapping symmetry. Furthermore, 
by requiring that H is everywhere non-vanishing inside the box, one gets the constraint 
X 2 < X 3 . In fact, a detailed analysis shows that to avoid a possible curvature singularity 
inside the box, we need only to consider the box satisfying the following ranges of coordinates 
and parameters: 

— 00<Xi<y<X2<X<X3< —X 2 , 

with X 4 obeying 

X 4 e (-CX), U (2^, CX)) , if X 3 < 0 ; 

X4e(2^,2^), ifX3>0. 

The special value X 3 = 0 can be recovered as a limit of either case. The two cases X 3 < 0 


(3.3) 

(3.4) 
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(b) 


Figure 1: The boxes for the cases (a) 0:3 = — 0:4 = |; and (b) 0:3 = |, 0:4 = 2. The other 

free parameters have been set as Xi = —2 and ^ The lighter-shaded strips represent 

the allowed ranges of x^. The dashed diagonal line represents asymptotic infinity. The solid 
curves are locations of potential curvature singularities, as are the crossed diamonds at the 
four points (x = 0, ?/ = ±cxo) and (x = ± 00 , ?/ = 0). Note that the shaded box in each case 
touches asymptotic inhnity at one corner while avoiding all the curvature singularities. 


and X 3 > 0 are illustrated in the plots in Fig. [T]^a) and (b) respectively, for specihc choices 
of parameters. 

When X 3 < 0, one can in fact further narrow the range of X 4 to the hnite interval 
—X 2 ) by using a combination of the inversion-swapping and scaling symmetries, which 

X\ 

maps a solution with any given X 4 = ra to one with X 4 = ^ while preserving the value of X 2 . 

The conditions that we have derived are also sufficient to ensure that the signature of 
the metric is positive- or negative-definite in the region of interest. Since we have required 
that X and Y have opposite signs, it follows that F defined in fl2.ip has the same sign as X. 
If k is positive, then the metric in fl2.ip can be seen to have a positive- or negative-definite 
signature. 

To summarise, we have the coordinate range X 2 < x < X 3 and xi < y < X 2 , which 
can be visualised as a box in an x-y plot. With the gauge choice fl3.2p . the solution is then 
characterised by five parameters: the position of the box determined by xi^ 3 , the fourth root 
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X 4 : of X, the line 


y = , 


(3.5) 


with a slope between —1 and 1, and an overall (positive) scale factor k. We have shown that 
these parameters have to satisfy (I3.3p and (13.dh . in order to ensure that the solution (12.ip 
has Euclidean signature and that it does not contain any curvature singularities. 


3.2. ALF and non-self-dual properties 

One key feature of the solution fl2.ip is its ALF property (with AF as a special case). 
One can quickly calculate that when x ^ y, 


1 

9 tt t - — , 

1 — 


(3.6) 


which is a hnite constant, a characteristic of an ALF metric. This selects a particular 
Killing vector ^ (up to constant multiplication) for the solution, which has a hnite norm 
and generates the Euclidean time how at inhnity. This property obviously indicates that, 
unlike say the Ricci-hat Plebahski-Demiahski solution, there is no symmetry between the 
coordinates r and 0. To show the solution is ALF, one can dehne the coordinates (r, 6 ) 
around the asymptotic region x = y = X 2 hj 


X2\/k{l - Z/2) 2^ X2\/k{l - Z/2) , .^9 

X = X2 -cos - , y = X2-\ -sm - , 

r 2 r 2 


(3.7) 


and (f, 0) via the simultaneous substitutions 


r —)■ \/l — f + 60, 0 ^ 


2 \fkx 2 


04(0:2 - Xi){x2 - X 3 ){X 2 - X4.) 


(3.8) 


For large r, with b chosen as 


b = 






{Xi - X2)(X2 - X3)(X2 - X4) 


the metric approaches fll.21) . with (r, 0) replaced by (r,0) and 


Vk[2i'X2{xiX2 + X2X3 + X2X4 — X1X3 — X1X4 — X3X4) — (1 + z^^) (0:10:30:4 — 0:2)] 

2^/T^^lF{x2 - Xi){x2 - 0:3)(0:2 - 0:4) 


(3.10) 
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or in an alternative form 


^ _ y/k[2u{aiX2 - - (1 - i^)‘^{ao - a^x^)] ^ 

2a4\/l - u‘^X 2 {x 2 - Xi){x 2 - X3){X2 - X 4 ) 

We refer to the Killing vectors associated to the newly dehned coordinates (f, 0) as the 
“natural Killing vectors” of our solution, which can be seen to generate natural notions of 
(Euclidean) time translation and rotation respectively at inhnity. In the limit n ^ oo, the 
solution should become ALE. Indeed, in the limits u = ±1, the solution reduces to the 
Ricci-flat Plebahski-Demiahski and triple-collinearly-centered Gibbons-Hawking solutions 
respectively, both of which are ALE. The solution can thus be thought of as a one-parameter 
family of ALE metrics interpolating between the latter two metrics, indexed by u. On the 
other hand, the solution becomes AF in the limit n = 0, and in particular it contains the 
AF instanton discovered in |T6] as a special case. These important limits will be discussed 
in detail in Sec. 4. 

The solution fl2.ip has a Riemann tensor that is in general neither self-dual nor anti¬ 
self-dual. One self-dual limit is taken by setting u = —1, which, as just mentioned, gives the 
triple-collinearly-centered Gibbons-Hawking solution. There is another self-dual limit in the 
case when X has two pairs of opposite roots, which gives the double-centered Taub-NUT 
solution. This limit will be discussed in detail in Appendix A.2. 


3.3. Rod structure 

The solution fl2.ip possesses two apparent Killing vectors ^ and It can thus be put 
in a canonical form, in so-called Weyl-Papapetrou coordinates given by 


_ yJ-XY _2{aQ\a2xy\a^x^'ip‘')R{x\y){axRazxy) (o 

f- (x-yr' 2 (x-yY ' * 

where z is determined up to a constant shift and a flip of sign. From the expression of p, 
one can see that the locations of the rods correspond to the roots (including possible inhnite 
roots) of X and Y. Turning points are then the meeting points of these rods: just like the 
Plebahski-Demiahski solution, the above solution contains three turning points in general. 

In the coordinates fl2.ip . the turning points are located at {x = X 2 ,y = Xi), {x = X 3 ,y = 
Xi) and {x = X 3 ,y = X 2 ) respectively, or equivalently, along the z-axis with 


2:1 = 


aA{xiX2 + x^Xa) 


Z2 = 


ai{xixj, + X2X4) 


Z 3 = 


a 4 {xiXi X2X3) 


(3.13) 
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respectively in Weyl-Papapetrou coordinates. The fonr rods have directions li = Ki/Ki, 
with 


= 

K 2 = 
K 3 = 

= 


—Qq + 2 uao + 2 uaiX 2 + 

2 ’ 
X 2 

z/^oo + 2na3xf + 2 i'aix\ — aix\ 


ry .2 

Xi 


—oo + 2z/ao + 2 i/aiX 3 + v a^x^ 

Xg 

z/^Oo + 2 va3x\ + 2 z/a 4 a;| — 040:2 


xr, 


,1 


,1 


Ki = 

K 2 = 

«3 = 


04 ( 3^2 - Xi){x 2 - X3){X2 - 3 : 4 ) 

2 Vkx 2 

04(0:1 - 3:2)(xl - 3:3)(xl - 3:4) 

2 \/kxi 

04(3:3 - 3:1)(3:3 - 3:2)(3:3 - 3:4) 


2 \fkx 3 


K4^ = Ki . 


(3.14) 


Here, the direction (a,/?) is dehned by a-^ + Alternative expressions for A'j)!] are 

Ai[l] = 04(3:13:33:4/3:2 + 2 u { xiX 3 + 3:43:4 + 3:33:4) — z/^o:^] , 

^2(1] = 04(3:4 + 2 u { xiX 2 + 3:43:3 + 3:43:4) - z/^3:23:33:473:1], 

^3(1] = 04(3:43:23:4/3:3 + 2 v { xiX 2 + 3:43:4 + 3:23:4) - y ^ x ^ , 

^'4(1] = 04(3:2 + 2 u{xiX 2 + X2X3 + X2X4) - z/^3:43:33:473:2]. ( 3 . 15 ) 


These are also usefnl, since onr analysis below will be based mainly on the factorised form 
of X in fl2.6p . involving its four roots. 

The above rod structure encodes much useful information about the solution. The rel¬ 
ative directions of the four rods, and the ratio of the lengths of two hnite rods are actually 
invariants of the solution. For example, the asymptotic NUT charge n is related to the 
directions of the two asymptotic rods by n = A(iF 4 (l] — 7^4 (Ij), where A is some proportion¬ 
ality constant. In the same spirit, one can define the NUT-charges carried by the individual 
turning points to be rz,* = A(A'j(l] — ATj+i)!]), so that one has 


_ \/k{Xi + VX2 Y{XiX2 - X3X4) 

1 ■ ^ 
23 :iV 1 - z/ 2(3:2 - 3:i)(3:2 - X3){x2 - 3:4) 

X 2 \fk{xi + VX 3 Y{x 2 X 4 - X 1 X 3 ) 

n 2 = - -j^=^= -, 

23:13:3 V 1 - v‘^[x2 - 3:1) (3:2 - 3:3) (3:2 - 3:4) 

_ \fk{x2 + VX3Y{x2X3 - X1X4) 

Tl^ — - / =- . 

2 XaVl - IX^{X2 - Xi){X2 - X3)(X2 - X4) 


(3.16) 


Note that the asymptotic NUT charge is then the sum of the individual NUT charges: 


n = Til+122 + 123. 

The rod structure will be extensively used in the subsequent study. In particular, many 
special cases of the general solution were identified by first studying the behaviour of the 
rod structure. In the rest of this subsection, we will briefly discuss the special case when two 
adjacent rods are joined up. 
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Two adjacent rods, say the i-th and i + 1-th rods, can be joined np by setting n* = 0. 
In this case, the turning point at which they meet is effectively eliminated, and one obtains 
a solution whose rod structure has only two turning points. For example, we can impose 

n3 = 0, (3.17) 

to join up the third and fourth rods, thus eliminating the third turning point from the rod 
structure. An obvious solution to this condition is 


XiX4^ = X2X3. (3.18) 

The resulting solution is the Kerr-NUT solution. In this limit, the box touches a solution 
curve of H at one of its corners. We have mentioned that the solution curves are locations of 
curvature singularities for general parameters. But in the present case, it can be explicitly 
checked that at this corner, a zero factor emerges from the numerator of the Kretschmann in¬ 
variant, which cancels the zero factor of H. This gives rise to a finite Kretschmann invariant. 
Another solution to the condition fl3.17l) is 


xs —)■ —X 2 , X 4 —>■ —xi , z/ — )■ 1, (3.19) 

which results in the double-centered Taub-NUT solution. This limit can also be understood 
as a box with one corner touching a solution curve of H, now corresponding to its factor 
ux + y. The details of how the Kerr-NUT and double-centered Taub-NUT solutions are 
recovered in these limits can be found in Appendix A.l and A.2, respectively. 

The joining-up of the first and second rods, or of the second and third rods, can be 
similarly obtained. We can see that the (finite) lower and upper bounds previously identified 
for the parameter x^ in fl3.4p correspond exactly to the joining up of different pairs of adjacent 
rods to obtain the Kerr-NUT solution. 


4. Various limits 

We have noted that the solution fl2.1l) is ALF, and that it contains two ALE limits: the 
Plebahski-Demiahski solution and the triple-collinearly-centered Gibbons-Hawking solution. 
In this section, we show in detail how these two limits are recovered. The AF limit of fl2.ip is 
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also discussed, and we show how the new AF gravitational instanton can be recovered from 
it. Finally, a new special case of fl2.ip is discussed: an ALF generalisation of the C-metric, 
which can also be called the NUT-charged C-metric. Note that all the limits discussed here 
have three turning points in their rod structure; limits with two or one turning points will 
be discussed in Appendix A. 


4 . 1 . Plebanski-Demianski solution 

Recall from fl3.6p that the norm of ^ at inhnity becomes unbounded when u —)• ±1. 
This means that the metric becomes ALE. Here, we will show that by directly setting 


1 / = 1 , 


(4.1) 


we recover the Plebahski-Demiahski metric from (12.ip . We hrst perform the Mobius trans¬ 
formation 

n-t-l « -U 1 

(4.2) 


p+l q+l 

X = - r , y = 


p — 1 ’ 


g - 1 ’ 


and redehne the parameters 


yem-t-n , , 

ao ,4 = TT - 7 T —^— , ai ,3 = T{m - n), 02 = - + 87 , k = -. (4.3) 

2 4 2 2 m — n 


Note that this redehnition of six parameters in terms of only four is consistent with the 
above-mentioned fact that two of the original parameters are redundant. Then after a linear 
transformation of the coordinates r and 0 via the simultaneous substitutions 


r —>■ 


(4m — 4n -|- 67 + e)r -|- (4m — 4n — 67 — e)0 


W' 


-)■ 


r 


m — n 


2^/m — n ’ 


(4.4) 


the metric is brought precisely to the Ricci-flat class of the (Euclideanlised) Plebahski- 
Demiahski metric [7] 


ds" = 


(p — g)^ L P 


-dp - 1 - 


P 


1 — p2g2 

Q 


1-pV^ 2 


dg" - 


(d 0 — g^dr)"^ 
-(dr -pM 0 )^ 


Q ^ 1 — p^g^ 

P = 7 -I- 2np — ep^ -I- 2mp^ -|- yp^, Q = 7 + 2nq — eq^ -|- 2mg^ -|- yg^. (4.5) 
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4-2. Triple-collinearly-centered Gibbons-Hawking solution 
The other ALE limit is taken by directly setting 


u = 


- 1 , 


(4.6) 


and one recovers the triple-collinearly-centered Gibbons-Hawking solution. To cast the 
resulting solution in a more familiar form, we hrst dehne the coordinates (r, 6 ) by 


r sin^ 


V-XY ^ 2(ao + a 2 xy + a^x'^y'^) + (x + y)(ai + a^xy) 

7 -7^’ rcos0 =---- 

[x-y)^ 2{x-yy 


(4.7) 


the parameters di, 2,3 in terms of the roots of X as 

a4{xiX2 + x^Xi) ayxiXs-i-X2X4) 


di = 


do — 


do — 


a^xiXi + X2Xs) 


(4.8) 


and ni^ 2,3 as 

_ k{xiX2 - X3X4) 

04 ( 0:1 - a: 3 )(a:i - 0 : 4 )( 0:2 - 0 : 3 )( 0:2 - 0 : 4 ) ’ 
_ k(xiX3 - X2X4) 

^ 04(0:1 - X2)(o:i - 0:4)(0:2 - 0:3)(0:3 - X4) ’ 

_ k{x2X3 - X1X4) 

04(0:1 - 0:2)(0:1 - 0:3)(0:2 - 0:4)(0:3 - 0:4) ■ 

After the simultaneous substitutions 


(4.9) 


r ^ (— 02 r + 4A:0)/\/4fc , 0^r/\/4fc, 


(4.10) 


the metric is brought to the form 

ds^ = V~^ (dr + A)^ + V (dr^ + r‘^d9‘^ + sin^ 9 dcp^), 

2 ni{r 


which is recognised as the triple-collinearly-centered Gibbons-Hawking solution. Note that 
there is a translational symmetry by adding an arbitrary constant to the right-hand sides 
of rcos9 in fl4.7p and <^ 1 , 2,3 in 114.81) . The charges ni^ 2,3 in fl4.9p always have the same sign 
in the ranges fl3.3p and fl3.4l) that we considered, as required by the absence of curvature 
singularities in the solution. 


—^cl 0 ^ rj = \/r2 - 1 - dy — 2dir cos 9 , (4.11) 

Ti 




i=l 


2 = 1 
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4-3. AF limit and the new AF instanton 


This limit is taken by setting the total NUT-charge flS.lip to zero, so now the metric 
CHi becomes AF. The solution to the condition n = 0 and the gauge choice X 2 = —1 can 
be written as follows: 

(1 — z/)^(ao — aA 

+ ^2 + ~ —- 2 - ■ 

Within this four-parameter AF subclass is the completely regular gravitational instan¬ 
ton recently discovered by the present authors [16] and christened the “new AF instanton”. 
It can be obtained as a special case by further imposing 


ii = £4 = ±£2 ± £3 ■ 


(4.13) 


These conditions ensure that the metric is defined on an underlying manifold with a globally 
well-dehned U{1) x 17(1) isometry, and is free of conical and orbifold singularities. A solution 
to these conditions can be written (for general X 2 ) in the following parametric forms: 


z/ = - 2^2 


Xi = 


^{l- 2 ^ + 2 e)x 2 

i- 2 e 


X 3 = 


(1 - 2 ^ + 2 e)x2 

4 e(i-o 


0:4 = 0, 


(4.14) 


or 


u = - 2 ^^ xi = 


4 e{i-ox 2 


X3 = 


(1 - 20^2 


0:4 = CX) , 


(4.15) 


1 - 2^ + 2^2 ’ «1 - 2f + 2f2) 

which are related to each other by the inversion-swapping symmetry. 

In the latter parameterisation with an appropriate gauge choice, the new AF instanton 
is then given by the metric fl2.ip with the parameter u and the structure function X defined 
as 


z/ = -2e^ X = {x + Ae- 4e')(a: + e - 2e' + 2e')(x - 1 + 20 . 


(4.16) 


The parameters ao ,...,4 are encoded in the polynomial X as its coefficients; in particular, 
we simply have 04 = 0 and a gauge-fixed 03 = 1. The solution is then determined by the 
parameters k and 0 This gives an alternative but simpler form of the new AF instanton. 
Now the awkward square-root terms in the latter are eliminated and its analysis based on 
this parameterisation will be much easier. To cast it in exactly the same form used in |16j . 
one needs to perform the Mobius transformation: 


0(1 - 2e + 20)^0: - (1 - 20)(1 - 4e + 20)] 
(1 - 2e + 2e)x + 1-20 


(4.17) 
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with a similar equation for ?/, followed by the simultaneous substitutions 


r 

k 




- 4^ + 8^^ - 12^3 + 16 ^^ - 8^^) ^ 

vT^^(i- 2 e)(i- 2 e + 2 e)' ’ 

4fc2 


^/r^^(i- 2 a(i- 2 e+ 2 a 

4fc^(l-2^(1-0^(1-20^ 

(i+ 2 a(i- 2 e+ 2 a" ’ 


and dually introduce the new parameters A and 7 defined as 


1 -2^2 ^ 1 - 4^ + 2^2 

1 - 2 ^ + 2^2 ’ ^ “ 1 - 2 ^ + 2^2 • 


(4.18) 


(4.19) 


NUT-charged C-metric 

The NUT-charged C-metric is obtained by requiring that the second rod is parallel to 
the fourth rod, which is achieved by setting 


04 = 0 or 0:4 = ±cxD , u = 0 . (4.20) 

The metric then becomes 

_ (^dr-F aoUd0)^ kH / dx^ dj/^ XU 2 ^ 

^ {x - y)HF {x-yY\X Y kF ^ 

H = y{-2ao - aiy + a 2 ,xyY , F = y^X - x^Y , 

X = oo + oix - 1 - a2x‘^ + a^x^, Y = qq + aiy + 02 ?/^ -|- a^y^. (4-21) 

We make the gauge choice 03 = 1 and use the roots of X = (x — xi)(x — X 2 )(x — X 3 ) as 
the parameters. Requiring H to be everywhere non-vanishing in the box gives the following 
ranges of the parameters and coordinates: 


— OO<Xi<|/<X2<X<X3<0, 


(4.22) 


which is consistent with fl3.3p and 03.41) from our general analysis. The reader is reminded 
that we actually have X 4 = —cxo, which provides the lower bound for xi and implies the 
upper bound for X 3 . 

The rod structure of this solution can be brought to a form very close to that of the 
C-metric by dehning 


T = T 


/ 




2 x2\/k 


(Xi - X 2 )(X 2 - X 3 ) 


(4.23) 
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( 4 n, 1 ) 




^(4n', 1 ) 


( 0 , 1 ) 


X = X2 z[ y = Xi 2;' X = X3 z'3 V = X2 


Figure 2: The rod structure of the NUT-charged C-metric, with asymptotic NUT charge n. 


In these coordinates, the corresponding rod directions are 

£i = ( 4 n, 1 ), £2 = ^(0,!), £3 = , 

K 2 H 3 

where n and K 2 are given by 


£4 = (0,1), 


(4.24) 


n = 


VkXiX3 


2 (X 1 - X2){X2 - X3) ’ 


Kr, = 


X2{xi - X3) 
Xi{x2 - X3) ’ 


(4.25) 


and n' and Kg are given by 


t .2 

X 2 


= T2 ^ > FT + 77 = 1 • 


X 


^3 ^2 ^3 

We note that the lengths of the two hnite rods are respectively 


(4.26) 


^21 ~ 


'/kX 2 

Xi - X2 


^32 ~ 


VkX 2 

X2 - X3 


(4.27) 


This rod structure is illustrated in Fig. 2. The first and third rods are the analogues of the 
horizon rods of the C-metric, while the second and fourth rods are the usual axis rods. 

Like the original solution, the metric fl4.2ip is ALF with an asymptotic NUT charge n. 
An ALE limit is obtained when n —)■ 00 , in which case we recover the C-metric. To take this 
limit, we need to set 

2 ^ 1,3 X 2 , (4.28) 

while hxing the other parameters. The box in this case thus becomes a point at {x = X 2 ,y = 
X 2 ). More specihcally, we hrst write the parameters in the following form: 


xi = 2:2(1 + (1 — c)e), X3 = ^2(1 — 2 ce), k = 


(4.29) 


and define new coordinates by performing the following substitutions: 


X ^ X2[l - c(l -F x)e ], y -)■ X2[l - c(l -F y)e ], 
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2xiX3Vk{l — c) 

T —)■ ^-r T , 




2 x 2 '/k{l — c) 


(r + 0). 


{XI-X2){X2-X3) ’ ' {XI-X2){X2-X3) 

After taking the limit e ^ 0, the metric becomes the familiar form of the C-metric: 


(4.30) 


= 


f 


H 

+ G{x)d4>‘^ —— G{y)dT‘^ , G(a;) = (1 — a;^)(l + ex). (4.31) 


(x- 2/)2 [G'(x) G{y) 

On the other hand, the AF limit of fl4.2ip is taken by requiring 


n = 0 


X3 = 0 or oo = 0. 


(4.32) 


In this limit, the hrst rod joins up with the second, or in other words, the hrst turning point 
vanishes. The resulting rod structure has two turning points, and turns out to be just the 
Schwarzschild solution. The mapping to this metric is a special case of the transformation 
described in Appendix A.3. 

The solution fl4.2ip admits several other, more subtle limits. In the limit Xi —)■ X 2 and 
/c —)■ 0 while keeping hnite, one can recover the self-dual Taub-NUT solution. In this 

case, the second and third turning points merge with each other and together they disappear 
from the rod structure. The details of how this limit is taken can be found in Appendix 
A.4. The self-dual Taub-NUT solution can also be obtained in the limit in which the second 
and third turning points are pushed to inhnity. These limits show that the asymptotic NUT 
charge can be thought of as being carried by the hrst turning point. 

Finally, in the limit xi = —00 or 03 = 0, the second rod shrinks down to zero length, 
and the hrst and second turning points merge with each other. In this case, one recovers 
the Kerr-NUT solution with NUT-charge parameter equal to the rotational parameter, i.e., 
flA.4p with n = a. 


5. Alternative interpretation as two touching Kerr-NUTs 

In this section, we show how the general solution fl2.ip can be regarded as a system 
consisting of two touching Kerr-NUTs. We also analyse in detail, the special case consisting 
of a Schwarzschild in superposition with a self-dual Taub-NUT. 

5.1. Kerr-NUT touching Kerr-NUT 

The general solution fl2.1l) has an alternative interpretation as a system consisting of 
two touching Kerr-NUTs: the south pole of one Kerr-NUT touches the north pole of the 
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other. This is consistent with our ISM construction. Recall that in our construction, one 
possible seed solution is the double Schwarzschild solution, with an inner axis separating 
the two black holes; the inner axis is then joined up with one of the black-hole horizons by 
appropriately choosing the corresponding BZ parameter. The inner axis disappears from 
the solution, and in this sense the resulting two Kerr-NUTs are touching each other, sharing 
one single point as their common pole. In what follows, we shall present the rod structure 
of the solution in a form that favours this interpretation, and show how one can remove one 
Kerr-NUT from the solution in certain limits. 

In the natural Killing coordinates, the rod structure is 

£i = (2n,l), £2 = -(l,f^i), 4 = —^'4=(-2n,l), (5.1) 

K\ Kii 


where the angular velocities hliji are given by 

\/l - u‘^XiX^{x 2 - Xi){x 2 - X3)(X2 - X 4 ) 
\/k(—Li -|- L2 + 

VT^^XiX3(x 2 - Xi)(x 2 - X3)(X2 - X 4 ) 


flj = 
hill = 


Ki = 

K\i = 


— L 2 + L 3 ) 

and the surface gravities ki n are given by 

\/l - u‘^X2X3{Xi - X2){xi - X3){xi - X4) 

\/fc(—Li -1- L 2 + L 3 ) 

yi- v‘^XiX2{x3 - Xi){x3 - X2){X3 - X4) 

■\/k{—Li — L 2 L 3 ) 

Here, Ti 2,3 are the three functions defined as 

Li = Xi{x2X3 - XiX 4 ){x 2 + 1^X3^, L2 = X2{XiX3 - X2X4){Xi -f ^3)^, 

L3 = X 3 {X 4 X 2 - X3X4){X4 ^2)^. 

We note that the lengths of the rods are 

y/k — ki>‘^X 2 {xi — X 4 ) _ \/k — kv‘^X 2 {x 3 — X 4 ) 

^21 = —7 -TT-) ^32 = 


(5.2) 


(5.3) 


(5.4) 


(Xl - X2)(X2 - X4) 


(X2 - X3)(X2 - X4) 


(5.5) 


This rod structure is illustrated in Fig. 3. For comparison, we recall that the rod structure 
of a single Kerr-NUT [12] has three rods, with directions (2n, 1), ^(1,11) and (—2n, 1) in 
suitable coordinates, where n is the NUT charge, and k and hi are the surface gravity and 
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(2n,l) (-2n,l) 

- m - m -•- 

X = X2 Zi y = Xi Z2 X = X3 Z3 y = X2 


Figure 3: The rod structure of the general solution, regarded as a touching double-Kerr-NUT 
system. 

angular velocity of the horizon respectively. It is then clear that the rod structure considered 
here consists of two Kerr-NUTs, whose horizons touch at the poles, kiji and fliji are the 
surface gravities and angular velocities of the two horizons respectively. 

One can remove a Kerr-NUT from the solution by zooming in to one of the Kerr- 
NUTs and at the same time sending the other pole of the second Kerr-NUT to inhnity. 
Alternatively, one can shrink the size of the second Kerr-NUT to zero. We will give a more 
detailed discussion of both of these possibilities in Appendix A.l. The explicit mapping of 
the hrst limit to the Kerr-NUT solution is given therein. Recall that in addition to these 
two possibilities, we have a third possibility to recover Kerr-NUT by joining up two rods as 
discussed in Sec. 3.3, corresponding to the lower and upper bounds of x^. 

In the rest of this subsection, we will briefly discuss a few special cases of this solution. 
First, consider the case 


xi = cxD, i.e., 04 = 0 III = 0, (5.6) 

in which the hrst Kerr-NUT becomes static]^ This can be interpreted as a conhguration in 
which a (non-self-dual) Taub-NUT touches a Kerr-NUT. On the other hand, the case 

X 3 = 0 , i.e., oo = 0 ^ On = 0 , (5.7) 

is essentially the same conhguration, with the locations of these two objects swapped. This 
is of course expected, since they are related by the inversion-swapping symmetry. 

The case with both Kerr-NUTs static: 


a:i = 00 , 0:3 = 0 , i.e., ao ,4 = 0 ^ = 0 , (5.8) 

can then be interpreted as a conhguration in which a Taub-NUT touches another Taub-NUT, 

which is of course a single larger Taub-NUT. We will discuss this limit in detail in Appendix 

^Note that there is a well-defined notion of being “static”, given by the requirement that the rod has a 
direction proportional to 
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A.3, and show how the familiar form of the Taub-NUT solution can be recovered. By further 
setting 1 / = 0 , one obtains the Schwarzschild solution, which can be viewed as resulting from 
a Schwarzschild touching another Schwarzschild. 


5.2. Schwarzschild in superposition with a self-dual Taub-NUT 

The case describing a Schwarzschild in superposition with a self-dual Taub-NUT is 
obtained by requiring that the second rod is static and that the hrst rod is parallel to the 
third rod. The solution to these conditions is 


2UX2X3, 

xi = —00 , X4 =-. 

0:2 -F X3 


(5.9) 


An equivalent solution can be obtained by applying the inversion-swapping symmetry, in 
which the third rod becomes static and the second rod becomes parallel to the fourth rod. 
To be concrete, we will focus on the former case in the rest of this subsection. 

Choosing the gauge 03 = —X 2 — X 3 , X 2 = —1, and dehning X 3 = p, the metric is then 
given by 

2 _(^dr + Gd0)^ kH fdx^ d|/2 2 ^ 

^ {x- y)HF ^ (x - vY 'V 

H = {ux + 2 /)[(l - /i)(ra - y)iix{ 2 iy - 1 ) - xy) + 4pV(l - u )], 

G = 2iy[{y\l - y) - - y{2p -l){x-yx- y)Y ], F = y‘^X- x‘^Y , 

X = {x + l){x — iJ,){x — yx + 2yi '), Y = {y + l){y — y){y — yy + 2yi') . (5.10) 


Recall that (x = 0 ,j/ = —cxo)isa curvature singularity, so one needs to impose y < 0. The 
requirement that H remains non-vanishing inside the box gives u > 0. Then the ranges of 
the parameters and coordinates are given by 


— l</i< 0 <i/<l, —00 < y < —1 < X < y. 


(5.11) 


These ranges are consistent with fl3.3p and fl3.4l) from our general analysis. The restriction 
on the range of the parameter u originates from its relation to X 4 in fl5.9p . 

The rod structure of this solution can be brought to a nicer form by dehning 


= \/l — R - 1 - 


Ayv{2v — l)^/k 
{I + y){l — y - 2 yp) 


2^^k 


{I + y){l - y - 2 yiy) 


(5.12) 
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(0,1) ;^( 1 , 0 ) ^( 0 , 1 ) (-4n,l) 

-•-^-•- - -•- 

X = X2 z'^ y = Xi z'^ X = Xj, z's y = X2 


Figure 4: The rod structure of a Schwarzschild in superposition with a self-dual Taub-NUT. 
In the parameterisation used here, one has xi = —oo, X 2 = —1 and x^ = fi. 


In these coordinates, the four rods are located at a; = —1, y = —oo, x = fi and y = —1 
respectively. The corresponding rod directions are 


£i = (0,l), £2 = -(1,0), 4 = -(0,1), £4 = (-4n,l). 


K. 


Kn 


(5.13) 


where 


Ko = 


\/r^ 




2 y/k 


K 3 = 


1 — fi + 2 ^ 


n = 


\/kv{l — fiiy)" 


I — fi — 2fiv ’ \/l — 1 / 2(1 + ^)(1 — fi — 2fiii) 

We note that the lengths of the two hnite rods are respectively 


(5.14) 


, ^ yfc(i - v^)(i - ji) , ^ yfc(i - 1 - 2v) 

l — — (1 + ^ — 2^1!^) 


(5.15) 


This rod structure is illustrated in Fig. 4. We recognise a Schwarzschild-like structure in 
the hrst three rods, while the third turning point can be regarded as an isolated self-dual 
Taub-NUT. The distance between the Schwarzschild horizon and the self-dual Taub-NUT is 
given by z'^ 2 - 1^ picture, the asymptotic NUT charge n is carried entirely by the self-dual 
Taub-NUT. 

Note that in the ranges of interest fIS.lip . the quantities K 3 , 2^21 z '^2 remain 

positive. The NUT-charge n is negative, and ranges from 0 down to —cx) when v ranges from 0 
to 1. We can thus interpret n as the asymptotic NUT-charge parameter: u = 0 corresponds to 
a zero value of the NUT charge and gives the Schwarzschild solution, while u = 1 corresponds 
to an inhnite (negative) value of the NUT charge and gives an ALE solution, k simply sets 
the scale of the solution and can thus be taken as the mass parameter, fi determines the 
length of the third rod, and so can be interpreted as the separation parameter: /i = — 1 
corresponds to the inhnite-separation limit z '^2 = gives either the self-dual Taub- 

NUT or the Schwarzschild solution; while fi = 0 corresponds to the zero-separation limit 
Z 32 = 0, in which the self-dual Taub-NUT merges with the Schwarzschild to give a non-self¬ 
dual Taub-NUT. We do not present the details here, apart from mentioning that some of 
these limits are special cases of those considered in Appendix A. 
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One could wonder what happens when — 1 < i/ < 0, and in particular when z/ = — 1. It 
turns out that these possibilities do not give rise to well-behaved solutions, since a curvature 
singularity appears inside the box. In the case v = —1, the above solution reduces to 
a triple-collinearly-centered Gibbons-Hawking solution with two nuts possessing opposite 
NUT charges. The latter is known to be singular, since it corresponds to at least one nut 
possessing a negative “magnetic mass”. 


6. Interpretation in Kaluza—Klein theory 

Given an ALF solution, there is a well-known procedure by which it can be turned 
into a solution of four-dimensional Kaluza-Klein theory. One hrst adds a flat time direction 
to the ALF solution to obtain a solution of hve-dimensional vacuum Einstein gravity with 
Lorentzian signature. Upon dimensional reduction along the compact direction ^ (or ^), 
one then obtains an asymptotically flat solution of Kaluza-Klein theory with Lorentzian 
signature. 

When this procedure is applied to the Taub-NUT solution, one obtains a magnetic 
monopole solution in Kaluza-Klein theory mi HE]- In this process, the NUT charge has 
turned into a magnetic charge associated with the Kaluza-Klein gauge held. When applied 
to the 77,-centered Taub-NUT solution, a system of n magnetic monopoles results. Each 
turning point in the rod structure of the original solution has turned into a source for a 
monopole. 

This procedure can also be applied to AF solutions with vanishing asymptotic NUT 
charge. The resulting Kaluza-Klein solution will then have zero total magnetic charge, 
although they would still in general describe a system of monopoles. For example, when this 
procedure is applied to the Kerr solution, one obtains a magnetic dipole na. 

It is clear then that the solution fl2.ll) can be turned into a Kaluza-Klein solution 
describing a system of three collinearly centered magnetic monopoles. Of the hve parameters 
in the original solution, three of them characterise the independent charges of the monopoles, 
and the remaining two their relative positions. The charges of these monopoles, as well as 
other properties of this system, can be read oh from the rod structure of the original solution. 
In the natural Killing coordinates (f,0), the four rods have directions U = Kj/kj, with 

ki = {2ni + 2n2 + 2n3,1), .A 2 = (-2ni -F 2n2 + 2n3,1), 

K3 = (-2ni - 2n22^3,1), .A4 = (-2ni - 2^2 - 2^3,1), (6.1) 
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where ni^ 2,3 ai'e the NUT charges of the turning points, given by fl3.16l) . Recall that the 
asymptotic NUT charge is n = ui + n 2 + ^, 3 . The surface gravities in these coordinates are 
given by 


3:2(a:i - X3){Xi - X4) . _ X2{X3 - xi){x3 - 3:4) 

Xl{x2 - X3){X2 - X4) ' ^ X3{X2 - Xi){x2 - X4) ' 


( 6 . 2 ) 


Now, ni, 2,3 can be identihed with the charges of the three monopoles. The identity 
= 1 indicates that the space-time described by the solution in Kaluza-Klein theory is 
asymptotically flat if (j) has standard periodicity 27r. Conical singularities then exist along 
the inner axes if \ki\ 7 ^ 1, for i = 2,3. For the ranges we consider, we actually always 
have 1 ^ 2 , 3 ! > 1, with equality holding only in certain degenerate limits. So in general, one 
cannot obtain a balanced triple-monopole system in Kaluza-Klein theory from fl2.ip . Conical 
singularities necessarily exist along the two axes joining the monopoles. This is in contrast to 
the triple-centered Taub-NUT solution, which describes a balanced triple-monopole system. 

Finally, we remark that in the AF limit n = 0, the three monopoles will carry a zero 
total charge. A special case of this system was previously considered in [TB] . 


7. Future work 

To summarise, we have presented a hve-parameter class of Ricci-flat solutions in four 
dimensions with Euclidean signature. This solution can be regarded as a generalisation of 
the Ricci-flat Plebahski-Demiahski solution with the inclusion of an asymptotic NUT charge; 
alternatively, it can be regarded as a system consisting of two touching Kerr-NUTs. Various 
properties and limits of the solution were studied. 

The reader should bear in mind that we have entirely focused on the solution as a 
Euclidean solution. The most prominent problem related to this solution is whether it admits 
a Lorentzian section and how to find it if it exists. Such a section would be very interesting, 
since it may have a (single)-black-hole interpretation, as the Plebahski-Demiahski solution 
does. We do not have a dehnite answer to this yet. On the one hand, the solution does 
contain various limits which do have Lorentzian sections, such as the Plebahski-Demiahski 
and Kerr-NUT solutions. On the other hand, the solution also contains the triple-collinearly- 
centered Gibbons-Hawking solution as a special case, which does not admit a Lorentzian 
section. This fact cannot rule out the existence of a Lorentzian section for our general 
solution, because a solution may have a Lorentzian section even though one of its subclasses 
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does not: as an example, the Kerr-NUT solution contains the Eguchi-Hanson solution as a 
special case and admits a Lorentzian section. We remark that if a Lorentzian section could 
be found, it is natural to interpret the lowest crossed diamond in Fig. [T]^b) as the black-hole 
curvature singularity, hidden behind the horizon sX y = Xi. 

It may also be interesting to study the algebraic properties of the solution. The 
Plebahski-Demiahski solution is algebraically special, and is so far the most general known 
Type-D solution in four dimensions. Based on the fact that our solution is a generalisation 
of the Ricci-flat Plebahski-Demiahski solution and that it possesses a very compact form, 
one may suspect that the solution may have some special algebraic properties. Of course, 
since a Lorentzian section is still lacking, one has to study this in the regime with Euclidean 
signature. 

Recall that the general Plebahski-Demiahski solution carries both electric and magnetic 
charges. A challenging problem is then to seek a generalisation of this solution, by adding 
electric and/or magnetic charges. It is even more challenging to seek an (A)dS generalisa¬ 
tion. We have mentioned that our solution emerges from the same ISM construction as the 
triple-collinearly-centered Taub-NUT solution. The latter solution has an (A)dS generalisa¬ 
tion, which can be obtained as a special case from the general construction carried out by 
Calderbank and Pedersen [TH]. These facts lead us to believe that it is worth more effort 
to consider the possibility of these generalisations in the future. One may even try to make 
guesses, by choosing an ansatz similar in form to our metric. 

The new AF instanton was identihed within this solution, as mentioned previously. 
Then a natural question to ask is whether the solution admits more new completely regular 
gravitational instantons. The analysis is straightforward with the aid of the rod structure 
presented in this paper, but the actual calculations are rather involved. More careful and 
exhaustive analysis is needed. 

Other possible future directions include the interpretation of the solution when embed¬ 
ded in string theory, in terms of a system of interacting D-branes. More detailed analysis of 
the solution in Kaluza-Klein theory is also worth pursuing. 
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A. Double-centered and single-centered limits 


In this appendix, we consider degenerate limits of the solution fl2.ip in which the final 
solution contains only two or one turning points. The resulting solutions include the Kerr- 
NUT, the double-centered Taub-NUT, the non-self-dual as well as the self-dual Taub-NUT 
solutions. 


A.l. Kerr-NUT solution 

It is known that a scaling limit [7] of the Plebahski-Demiahski solution results in the 
Kerr-NUT solution (as a subclass of the Carter-Plebahski solution [Sill]). This limit corre¬ 
sponds to sending the first or third turning point to infinity while zooming in to the region 
around the remaining two turning points. This limit is still present in the general solution 
fl2.ip . although in a more subtle way. To recover it, one first defines the parameters 

xi=X 2 - 2 e, xs =-X2 - cie , X4 =-X2 - C2e , z/= 1 -F C3e/x2 , k = l/e^, (A.l) 


and the coordinates (r, 6) by 
Arx2 


X = x^ + e 


■X2 


+ 


(ci -t- C3)(2 -|- Cl -f C3) 


/(2 -|- Cl -f- C2 -|- 2C3) 2 -|- Cl -|- C2 -|- 2 c3 

After making the substitutions 

a: 2 V^(i^ 2 [l] + i^4[l]) 


, y = X 2 — e{l + cos 6). 

(A.2) 


r 


-)■ 


e(2 -|- Cl -1- C2 -|- 2C3) 
X2 

2X2^/^ 


r 


a^ixi - X2){X2 - X3){X2 - X4,) 


(A.3) 


ai{xi - X2){X2 - X3){X2 - Xi) 

and taking the limit e —)■ 0, one can show that the first turning point is pushed to infinity. 
The metric in this limit becomes 

A ri ^ ^ 

ds^ = —[dr -I- [2ncos9 -|- asin^6') d0]^ -I-—— [ad^jJ — (r^ — — a^) d0]^ 


T, = — {n — a cos 9)\ A = — 2mr — + n^, 


(A.4) 


with parameters 


m = 


\//[(ci -f + (c2 -f 03)^ -|- 2ci -|- 2 c2 -|- 4C3] 
V-64a:i(2 -|- Cl -f- C 2 -|- 2 C 3 ) 


















n = 


V^[(ci + C3)(C2 + C3) + 2 + Cl + C2 + 2 C 3 ] 
1 ^/—T 6 x|'( 2 ''+~c 7 ^i-~C 2 ~+'^c^ 


a = 


/(2 + Cl + C2 + 2C3) 


16 x 9 


(A.5) 


This is the familiar form of the Kerr-NUT solution. 

However, there is a more natural limit in which one can recover the Kerr-NUT solution. 
As mentioned above, one can join up two adjacent rods to eliminate the turning point 
between them, instead of sending it to inhnity. Without loss of generality, here we consider 
the limit fl3.18p . in which the third rod is joined up with the fourth. We first define k in 
terms of I, and choose 04 appropriately: 




2 l{xi -X2){Xi -X3) 
xi{x2 - X3) 


04 = 


4/ 

{Xi - X4){X2 - X3) ■ 


(A. 6 ) 


We then bring the solution to its natural Killing coordinates by doing the substitutions 


T \/l — — b(j )), 0 —)■ —(f )/\/l — 1/2 , 


(A.7) 


with b given by 


^ _ 2/[(l - v‘^){x\ + x\) + 2v{{xi -t- 0:4)(0:2 + X3) + 2xia;4)] 

(1 - v‘^){xi - Xi){x2 - X3) 


(A. 8 ) 


By introducing the parameters m, a and n 

_ l[{l + v‘^){x2 + X3) + 2 v{xi + Xa)] _ 1 {xi-xa) 
^ (1 - u‘^){x2 - X3) ’ “ a;2-a;3 

_ 1[{1 + u‘^){xi + Xa) + 2u{x2 + X3)] 

(1 - U‘^){X2 - X3) 


(A.9) 


and the coordinates r and 6 


r 

cos 6 


l{x 2 + X 3 ){x + y) - 2 {xy + X 2 X 3 ) 1{1 -h i^‘^){x 2 + X 3 ) + 2 i/{xi + xa) 

(x2 -X3)(x-y) (1 - iy^)(x2 - X3) 

{xi -h xa){x + y) — 2 {xy + xiXa) 

{xi-XA){x-y) 


the metric is brought to the Kerr-NUT solution in the familiar form flA.4H . The corresponding 
joining-up limit of the Plebahski-Demiahski solution results in the double-centered Gibbons- 
Hawking solution, a subclass of the Kerr-NUT solution with inhnite asymptotic NUT charge. 
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A third way to eliminate a turning point from the rod structure is to set the length of 
say the second rod to zero, while keeping the third rod finite. In the touching double-Kerr- 
NUT picture, this corresponds to shrinking the size of the first Kerr-NUT to zero. This limit 
is achieved by taking 

03,4 = 0 or ao,i = 0 . (A.11) 

We will not present the detailed coordinate transformations here. 


A.2. Double-centered Tauh-NUT solution 

To recover the double-centered Taub-NUT solution, it is convenient to use the natural 
Killing coordinates (r,0) defined in fl3.8p . We define 


u = l — cie/l, k = 2cil/e, X 3 =—X 2{1 — e), X 4 =—xi{l — C 2 e), (A.12) 


and take the limit e —)■ 0. We then introduce the new parameters (a,ni,n 2 ) by 

_ 2ciXiX2 _ 1{1 + C 2 ){xi + X 2 ) _ 1(1 - C 2 ){xi - X 2 ) 

xj — xl^ ^ A{xi — X2) ’ A{xi -h X2) 

and the coordinates (r, 9) by 


r sin 9 
r cos 9 


2c\\/—{p? — x\){x‘^ — a;|)(|/2 — x\){y‘^ — x^) 

{xj - xl){x - 
2 ci{xy — x\){xy — xD 
{xl-xl){x-yf 


(A.13) 


(A. 14) 


After dropping the tildes on the natural Killing coordinates, and taking 0 —cj), the metric 

is brought to the familiar form: 


ds^ = V ^ (dr -I- A)^ -|- V (dr^ -|- r^d^^ -|- sin^ 9 dcfp '), K = 1 -|- , 

tr 

A = —— -^cl0, Ti = \/+ di^ — 2 dir cos 9 , di 2 = To . (A. 15) 

^ r,: 


The parameter ranges that we are interested in, fl3.3p and (13.41) . imply that —1 < C 2 < 1, 
which in turn ensures that Ui 2 are non-negative, consistent with the requirement that the 
“magnetic masses” are non-negative. 
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A.3. Non-self-dual Taub-NUT solution 
As mentioned in Sec. 5.1, the case 


flo,4 — 0 , 


(A.16) 


describes a (non-self-dual) Taub-NUT touching another Taub-NUT, resulting in a larger 
Taub-NUT. To map this solution to the familiar form of the Taub-NUT solution, we choose 
the gauge 03 = 1 and write X = {x — X2)x{x — xf). Note that one has xi = —00 and X 3 = 0. 
We then define 


r = 


'/k{u'^x — y) 


y/l - z/ 2 (x - y) ’ 
and make the substitutions 

r —)■ \/l — v'^T -|- 


cos 9 = 


2{xy + X2X4) — {x2 + X4){x + y) 


(x2 - Xi){x - y) 


2i/'/k{x2 + X4) 


-)■ 


2^/k 


X2 — X4 X2 — X4 

The metric then becomes the Taub-NUT solution, in the form flA.4p with 


(A.17) 


(A.18) 


\/k{l -|- z/^) uy/k 

2y/\ - z/2 ’ 


a = 0. 


(A.19) 


It is interesting to note that, in this limit, ^ 2,4 make no appearance in the physical 
parameters and become redundant. One is thus free to choose any value for 01 , 2,3 without 
changing the physical space, as long as the solution does not become degenerate. The non¬ 
trivial parameters are k and v, which can be interpreted as the mass and NUT charge 
parameter respectively. It can be seen that the special case with oq = 04 = = 0 is the 

Schwarzschild solution. 

This is also an example of joining up two rods: the second rod is joined up with the 
third. If instead one insists that —00 < X 4 < X 2 , then one has to close the box by x = X 4 
(instead of x = xi) as its second rod. In this case, it is the hrst rod that joins up with the 
second. The hnal solution is of course still the non-self-dual Taub-NUT solution. 


A. 4 . Self-dual Taub-NUT solution 

In this subsection, we show how the self-dual Taub-NUT solution is recovered from the 
NUT-charged C-metric fl4.2ip . We first dehne the parameters and coordinates 

xi = X 2 — e , k = le‘^, X = X2 + we , y = X2 — ze , (A.20) 
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and make the substitutions 


r ^ r + 


y/kxiXs 


{Xi - X2){X2 - X3) 


-)■ 


2\/kx2 


{Xi - X2)ix2 - X3) 


(A.21) 


After taking the limit e —)■ 0, the metric then becomes 

ds^ = V~^ (dr + 2n cos 9)'^ + V (dr^ + r'^dO^ + sin^ 9 d0^), V = 1 -\ -, (A.22) 

r 

if one identihes 


r = 


Vlx 2 {z — w — 1) 


w + z 


, cos 9 = 


w'^ + w + z"^ — z 
w'^ + w + z — z'^ ’’ 


n = 


VIX2X3 
2 {x 2 - X3) 


(A.23) 


This is the standard form of the self-dual Taub-NUT solution. 
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